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Abstract 

We investigate the implications of i?-Parity violation (RPV) for quark flavour viola- 
tion both by constraining the sneutrino masses to be positive and by studying the 
processes 6 — 57 and — mixing. In the latter there are two major contribu- 
tions, one from "direct" one loop diagrams involving RPV couplings, and one from 
the "indirect" contributions generated by the renormalisation group. We compare 
the effects and discuss the implications of our results. 



1 Introduction 



One of the most promising candidates for physics beyond the so-called Standard 
Model (SM) is that of supersymmetry (SUSY) [|l|. SUSY has the highly attractive 
properties of giving a natural explanation to the hierarchy problem of how it is possi- 
ble to have a low energy theory containing light scalars (the Higgs) when the ultimate 
theory must include states with masses of order the Planck mass. This allows the 
ultimate hope of constructing a theory of gauge unification and fermion (and super- 
partner) masses defined at a scale near the Planck mass yet whose structure and 
parameters may be deduced from physics at accessible scales. 

In this paper we shall be concerned with the implications of a particular possible 
feature of SUSY, namely that of i?-Parity violation (RP V) Q . i?-Parity is a "Z^ 
symmetry of both the SM and its minimal SUSY extension, the MSSM, under which 
all of the SM particles have charge 0, while all their SUSY partners have charge 1. Its 
implications include the stability of the lightest supersymmetric particle (LSP), and 
hence the typical SUSY collider signatures of missing Eq^ and the existence of a source 
of dark matter. Its violation changes both the implied cosmology and the expected 
collider signatures, allowing such effects as LSPs decaying inside the detector and 
leptoquarks 0. In addition to these, further constraints on RPV can be derived by 
considering experimental limits on rare decays |^, § , and by the demands of proton 
stability. In practice it is usual to evade the problems of proton decay by considering 
either baryon number or lepton number violation but not both simultaneously. 

i?-Parity is violated by the superpotential and soft potential 

W = h'l^Q,H2U, + h'}^Q,H,d, + h'r^L.H.e, 

+—XijkLiLjek + X'ijkLiQjdk + —X'-jj^Uidjdk 
+fi4^HiH2 + fiiLiH2 
Koft = VijQiH2Uj + r]fjQiHidj + rj^jLiHiej + h.c. 

1 



~^'2~'ijkLiLjek + C^jj^LiQjdj^ + —C'^jf^Uidjdk + h.c. 

+D^H,H2 + DiLiH2 + h.c. (1.1) 



From the point of view of deriving constraints on the i?-Parity violating couphngs in 
the model, the most extensively studied couplings are the dimensionless couplings A, 
A', and A", which directly generate many effects which can be experimentally limited. 
The extra soft terms by definition mostly couple only heavy SUSY particles and hence 
are relevant mostly because of their impact on the Renormalisation Group Equations 
(RGEs), although they can have significant effects on the neutrino- neutralino and 
Higgs boson-sneutrino sectors |^, p|, p!0 . 



In SUSY models, flavour- changing effects may be caused by the existence of off- 
diagonal terms in the sfermion mass matrices in the basis in which the fermion masses 
are diagonal, in which there are no tree level SM flavour- changing neutral currents 
(FCNC). Such flavour- violating soft masses can be generated either from the high 
energy theory such a GUT directly, or else through the RGEs by couplings which 
violate flavour symmetries, such as Yukawa couplings mediated by the CKM matrix 
and both possibilities have been studied extensively |Tl|, |12[. However, the inclusion 



of the RPV couplings in the RGEs allows many extra violations of quark and lepton 
flavour. 

In a previous analysis ||10|, we presented the renormalisation group equations 



(RGEs) for the couplings of the full i?-Parity violating sector of the model, and 
investigated the implications of typical scenarios at the GUT scale for the generation 
of neutrino masses and the decay fi 67. Contributions can conveniently be split 
into two categories which we term "direct" (where the flavour violation occurs directly 
through i?-Parity violating vertices in the diagrams) and "indirect" (where i?-Parity 
violation induces flavour violation through the RGEs on the soft masses). We found 
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that the indirect effects were typically large, often orders of magnitude larger than 
the direct ones. However, the extremely complicated dependence on the spectrum 
and the existence of many possible cancellations in the amplitude render the process 
of deriving bounds impossible except on an order of magnitude basis. 

The intention of this paper is to extend the work of reference [l^ to consider 



quark flavour violation (QFV) effects (for recent work on the subject see ref. [13|) . 
Following this introduction, we give a short discussion on a new bound from requiring 
the sneutrino masses to be above their experimental limits in section 2, while section 
3 is devoted to the effects of RPV on 6 ^ s'j, and section 4 to the — mixing 
term Amx- Many useful formulae for our RGBs and definitions are relegated to the 
appendices, while section 5 contains our conclusions. 

2 Bounds from Sneutrino Masses 



We now comment on a new simple bound on certain of the Yukawa couplings, which 
is associated with our assumptions of unification of couplings at a high scale. If we 
add a new Yukawa coupling associated with the sneutrinos (here A'), the mass of the 
sneutrino is driven down by its effect on the RGEs, and so increasing the value of 
the coupling will ultimately drive the mass of the sneutrino below the experimental 
limit. The low energy values of the lepton masses and the corresponding soft masses 
are given by 

ml = ml + 0.51M^/^-Y.X%,{MGUTmml + A9M^/^-1.5M,/2Ao-12Al] 

jk 

ml = ml + ^Mlcos2l3 (2.1) 

where we have assumed universal masses at the unification scale and solved the RGEs 
numerically. The experimental limit on sneutrinos is 37 GeV [|14| (assuming that the 
sneutrinos are not degenerate), so that we conclude that the Yukawa couplings A' 
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must be bounded by 

,2 m2 + 0.51M2/2 + |M|cos2/3-(37GeV)2 

1. A .,,iMauT) < i3^2 + 49^^^_i.5^^/^^^_i2Ag (2-2) 

Given that masses for squarks and sleptons as low as 100 GeV are barely tenable in a 
unified framework with present collider limits, this constraint is the tightest available 
on certain A'. For example, if we set Aq ~ 0, M1/2 = itlq = 200 GeV and tan/5 = 10 
then this gives a bound on all of the A' of 0.15, while if two of the A' are equal then 
they must be less than 0.11, bounds which become only very slightly weaker with 
increasing soft masses, while the corresponding numbers with mo and M1/2 both 100 
GeV are 0.12 and 0.09. Note that these couplings are at the GUT scale, and the 
electroweak scale values are then around three times larger, while there is an error of 
at least 10 to 20% from the uncertainty in the value of the strong coupling. 

Although we have neglected the effect of the tau Yukawa coupling, this in fact only 
makes the bound for the third generation rather tighter by reducing the sneutrino 
mass for that generation still further. We can also derive bounds on the A and A" 
couplings from similar arguments, but the results are not very restrictive. 

3 6 ^ S7 

3.1 Contributions 

The rare decay process 6 ^ 57 has a branching rate which can be deduced from the 



decay B K*'j which has been measured |T5[ to obtain a value 

B(6 ^ 57) = (2.32 ± 0.57 ± 0.35) x 10"^ (3.1) 
consistent with the standard model result, and hence a 95% confidence level limit of 

1.0 X 10"'* < B{b 57) < 4.2 X 10"* (3.2) 
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This gives a very strong test of new physics such as supersymmetry, since SUSY 
models with hght spectra can give large contributions to this process through charged 
Higgs and chargino diagrams 

The total branching ratio for the process b —>■ sj can be written (in units of the 
BR for the semileptonic b decay) as: 

B(6^g7) _ Sna / - . ~ 
B{b cev) Gj,\Kcb\^I[z) ^ ' 

where Gp is the Fermi constant, z = mc/mb, I{z) = 1 — Sz"^ + 8z^ — 24z'^ log(2;) is the 
phase space factor, Kij will be the different CKM matrix elements, Alr, Am are the 
total amplitudes to the LR and RL transitions respectively, and 

contains NLO effects {n{z) being the NLO correction to the semileptonic decay). Here 

A, = r^l6/23^7 + 8 (^14/23 _ ^16/23)^9 _^ ^^0 (3 5) 

o 

with i = LR, RL, rj = as{Mw) / Ois{mh) and the different terms are as follows: ^lrrl 
are the coefficients of the effective operators for the 657 and bsg interactions 

CIr,rl = ^rnb{sa^'^Pn,Lb)F,, 



CIr,rl = ^M^^'^'PrMGI^T^^ (3.6) 



AP^ = —awK^gKtb/M^ is from the coefficient of an operator 

= [sPLc][cPLb] , (3.7) 

and C stands for the leading logarithmic QCD corrections (for a complete list of 
references see [|T7l). 



The LR amplitudes can be divided into an i?-Parity conserving part plus an RPV 
one; the former has been calculated in ref. |T^, and corresponds mainly to the contri- 



butions coming from the SM diagram plus those with top quark and charged Higgs, 



and stops/scharms and charginos running in the loop, plus smaller contributions from 
loops with neutralinos or gluinos and d-type squarks, which are generated due to QFV 
explicitly through the CKM matrices. Their expressions are 



Al,9 



Al>9 



2 *^ *V2,M2,^^'nM2 



w , 




1,9 



m 



(3.8) 
(3.9) 




y/2Mw sin 13 













j 







(3.10) 



t—weak 



'^i—diag 



where our notation is as in ref . |jTO[ . T, B are the orthogonal matrices that diagonalise 
the stop and sbottom mass matrices respectively through TMj 
where Tn = T22 = cos T12 = —T21 = sin^^, so that we may write the mass eigen- 
states > and > as cosOtlii > +sin6't|t/j > and — sm9t\iL > +cos9t\iR > 
respectively, and similarly for other flavours. For the first and second generations 
with small left-right mixing we shall take cos6'j = 1, so that for example c/-^"* = d]^ 
and c/-^"* = c^. For simplicity, we shall use the notation §1 rather than d'2 \ except 
where confusion might arise. The different functions /^*^ are defined in Appendix C. 

We generate new contributions to both the LR and RL amplitudes from i?-parity 
violating couplings (A' and A") directly and also indirectly from the induced QFV 
soft terms. Therefore we can write 



^LR — ^LR + ^LR 



-Am _ -Am n rAm- 
^ -l_ A ^ _|_ A S 



^LR 



(3.11) 
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The direct amplitude is given by 



477 



12 



sin^ cos^ 6'rt 



+ 



r{l) 



m 



r(2) 



(3.12) 



where A' are the different RPV couphngs in the fermion mass eigenstate basis which, 
in this case, is related to the weak one by: 



A ijk A imkKm.i 



(3.13) 



where a sum over m is understood. Note that this still leaves the possibility of 



generating effects even with only one non-zero RPV coupling in the weak basis [[19 
Here Xji = ml. /ml.. 

The other amplitudes are: 



A ^ 



^LR 



^LR 



2 - 2 



n=l ^tr. 



m 



n 



n2 



ci j=l 



X 



2 



ml 



v/2Mvt/sin/?, 



(3.14) 



^2 cos 13 Mw 
2 Am^^^ 4 
2ai4^(5d E 712 — ^TTF E 



n=i - m-s, 



swQdN'ji 

- (swQdN'^, - ^(1/2 + Qds'^wWj,^ 



^ il/2 + Q,sl,)N'.^' 



Cw 

1 

Cw 

mb 



F2(X 



2 

mf 

On 



m\ 



swQdN'^i - -^QdN'j2\ K2 



2Mw cos (3 



mb 



2 

mf 

On 



m\ 



(3.15) 



2 Am2 - 
-2«,Q,C(i?) — 



2 

1 mi 



X 



B. 



F2(x'' 



gnJ 



F2{xl 



nl 



2 

On 



^n2 



F4(X^ 



Fiixl 



9^' 



ml 



2 

mi 

On 



m\ 



(3.16) 



2 Am2 - 
E 



-1 



m 



-B*aiC{G) 



Fiixl) F,{x 



31 



m~ 



m 



+S:i(2C(i?)-C(G)) 



m? 



m 
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-B;,^(2C(fi) - C(G)) - ^ 1 J , (3.17) 



where is the strong gauge couphng constant. Also, 



7I 2 ' 7I 2 ' q1 2 

^•''^ ^ rrd ' ^ ^ ^ ^i" ^ ^ ^' ^gn^ ^^-^ (3.19) 

^ii 

and the QCD factors are C{R) = 4/3, C{G) = 3. 

Now we turn to the RPV contributions to Arl. These are given by: 

= a'rl + A';, + i^r^" + , (3.20) 




where both A' and Xji are defined after eq. ( p.l2|) . 



A 



RL 



-2E 



3 \" \" 



An 



sin^ e. 



sin^ e, 



.f(5) 
J -t.a 



■{1) 



a) 



cos^ 9(1 . 



(2) 



f(5) 
•'7.9 



m: 



Mi 



3 



m 



■''7,9 



(1) 



m 



(1) 



cos^ 6. 



•I 7,0 



QdA 



RL 



Am? 



S2bn 



-2awQd / . 2 2 

^1 - i=i 



E 



(3.22) 



Cw 



X fi. 



n2 



^2(0;^ 



m 



S2 



- {swQdN';,--^QdN;i 



Cw 
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IgAm-g 



X 



X 



swQdN!jl - —(1/2 + Qdsl)N!^l ) Ki + 



rrib 



j2) 



2 

On 



-2asQ,C{R) J2 



Am? 



^n2 



m 



-a. 



E 



Am? 



m 



S2 



m 



S2 



m 



m;, 



2 



bn 



m 



S2 



m 



S2 



mfo \ mf 



m 



S2 



(3.23) 




(3.24) 



(3.25) 



3.2 Analytical Discussion 



In general the results for the implications of the MSSM for b sj are well known [[16 1. 
The SM contribution to the branching rate is large and in good agreement with the 
data at around 3 x 10""^. The MSSM charged Higgs contribution is always of the same 
sign as that from the W and can easily push the total over the experimental limit, 
hence giving quite a tight constraint on the charged Higgs mass or equivalently m^. In 
practice, however, the situation changes when the full spectrum is considered, because 
the SUSY partners can give contributions to the amplitude of both signs with the 
dominant contribution being that of the chargino which can be of opposite sign to that 



of the SM and charged Higgs |]16|. There are thus two separate scenarios according 
to the sign of /i. For /i4 < 0, the structure of the chargino mass matrix and mixings 
makes its contribution comparable in size to the other two with the opposite sign 
giving rise to a strong cancellation, and hence the constraints on SUSY parameters 
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are rather weak. For /i4 > the chargino amphtude is not so big and moreover its sign 
is not always opposed to that of SM and charged Higgs, and therefore cannot cancel 
them off completely so we must have a fairly heavy charged Higgs, which given our 
unification and radiative electroweak symmetry breaking assumptions implies fairly 
large soft masses. 

If we now turn to the implications of i?-parity violating couplings, there are three 
distinct possibilities. The first of these is where the product X'i2j^'i3j is non-zero. 
Here the main effects will be to change Am, since mixing is generated in left handed 
superfields only. This will give an additional term which is added directly to the 
amplitude from the usual MSSM terms, but may be of opposite sign and hence less 
restricted. By contrast, the cases of non-zero A'jj^A'ijs or A-'^-gA'/js generate contribu- 
tions to AjiL (negligible in the MSSM), which cannot therefore have interference with 
the MSSM effects but which are only added to them in quadrature. 

Beginning with the direct i?-Parity violating contributions, the SM contribution 



to Alr is — 6.7x 10 ®GeV ^, with an experimental limit on y {A-Ir}'^ + (^rl)^ derived 
from equation (U) of 4.1 to 8.5xlO"^GeV"l Since the QCD corrections typically 



give a contribution to Ai of around 0.65A], we shall consider what values of the 
couplings give a contribution of -2 x 10-^GeV"^ to Al^ or 7 x lO-^GeV"^ to A^^ 
(recall that the sign of the former matters, since it is being added to a negative SM 
contribution, while the latter is added in quadrature). It is then straightforward to 
derive the bounds 



A'i2jA'i3i(Afz) 

|A'.,2A',,3(Mz)| 

\K2jKdMz)\ 



~ 0.09 



0.035 



~ 0.16 



( Y 



)■- 


( lOOGeV^ 


2 








2 




2 




< lOOGeV^ 






V ^4 J 





(3.26) 

(3.27) 
(3.28) 



VlOOGeV/ 

where we have neglected the distinction between A' and A' in deriving the second of 
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these. 

The indirect contributions are far more comphcated, but it is instructive to con- 
sider their size in the case where we include only the diagrams with a helicity flip on 
the gaugino line (and hence an overall factor of the gaugino mass), set all sparticle 
masses approximately degenerate at m and ignore the various mixing factors. This 
then gives an approximate contribution for each chargino of 

Alr ~A.2,A.3,(Mgut) ^^^^ [-^^) ^^-^^^ 

where we have used equation ( p.l|) and have evaluated the functions with all masses 
equal. Since the product of two of these couplings at the GUT scale is an order of mag- 



nitude smaller than at Mz for X'i2jX'i3j we expect a significant chargino contri- 
bution relative to the direct contributions for large tan (3 unless the gaugino-higgsino 
mixing is small. A simple way of understanding how the chargino contribution must 
be relevant is to consider this in terms of results from the MSSM. The chargino con- 
tribution here is driven by the CKM matrix element K23 — 0.04, while the i?-Parity 
Violating contributions are driven by a similar mixing Am^/m?, which is roughly 
equal to the product of X'i2jX'i3j{MGUT) and hence will be similarly important when 
this product of couplings at the GUT scale is around 10~^. 

The gluino contribution is simpler, and neglecting mixing between sbottom squarks 
and the mass difference between squarks and gluinos we find a bound on each of the 
pairs of couplings 



A'.,2A'.,3(Mz) ^ 0-006 (y^^)' (3.31) 



K.,Ks,iMz) ^ 0.006 (y^)' (3.32) 

where the bound is in fact on the sum of these indirect terms plus the direct con- 
tributions quoted above. Unfortunately, we might expect the suppression from the 
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mixing in the sbottom sector to weaken these by an order of magnitude or so. We 
can conclude however that it is hkely that the RL contributions are dominated by 
the gluino, especially in the case of non-zero A^^^-^ijs; while for X'i2j^'i3j the chargino 
is likely to dominate, results which are confirmed by our numerical studies. 

3.3 Numerical Results 

Given the discussion above, we are now in a position to discuss the effects of set- 
ting these products non-zero with a realistic spectrum from universality at the GUT 
scale. The products of couplings which we are considering have not been constrained 
before, but we find X'i2j'^'aj < 0.06(100GeV/m)^ using the bounds on each coupling 
independently from decays |19|. Other such bounds can always be evaded by 



selecting the indices in the products appropriately; for example the tight bounds on 
X'ljf^ from neutron anti- neutron oscillation ^ only apply if i = 1. 

In Figure Q we show a simple example of mixing involving the product X'i2jX'i3j, 
where we choose both couplings 0.05 at the GUT scale, hence making the product 
around 0.023 at low energy, and with other parameters tan/5 = 10, /i4 < 0, mo = 100 
GeV, Aq = 0, and M1/2 varying. We see that in fact here the gluino contribution 
is dominated by the direct one, and this in turn is dominated by the chargino con- 
tribution. However, even at very low values of the mass spectrum (for M1/2 ~ 110 
GeV the chargino mass is actually below the experimental limit) and even though 
we have taken the product of couplings to be close to the limit of equation (p^), the 
contributions are still too small to greatly affect the total. Hence we conclude that 
bounds on this product are really not usefully found from 6 — > 57 decays. 

Similarly, for X'ij2X'ij3 we expect that the product at the GUT scale will be ~ 0.01 
from equation (|2.2|) , which in turn leads to the product at the weak scale being around 
0.1, which can only give a useful bound if the spectrum is extremely light, and in 



< 
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practice it again appears that constructing a spectrum so light as to have interesting 
consequences for 6 — > 57 through i?-Parity Violation requires one or other of the 
masses to become lighter than its experimental bound. 

For the case of A^j^A^jg we do not have such tight bounds from the sneutrino mass 
and can increase the coupling to quite a large value. In Figure we show a plot with 
the same input parameters as in Figure |I|, but with ^ij2Kj3i^GUT) = 0.01. Here 
again the effects on 6 — 57 are only really beginning to be significant for a very 
light spectrum. However, it is interesting that the contribution is very much larger 
than might be expected from equation (|3.28 ), since the gluino term is overwhelmingly 
dominant. 

In conclusion, while the i?-Parity violating couplings give a contribution to 6 — 57 
which is enhanced typically by an order of magnitude when the indirect contributions 
are included, and although this in principle gives a tighter bound than others in the 
literature, in fact these constraints are very weak in the context of universal mass 
spectrum at the GUT scale. This is because, for example, it is virtually impossible 
to arrange that squarks should have masses as light as 100 GeV without violating 
one or another of the experimental limits. In the case of non-zero A' couplings the 
bounds from 6 — 57 are weaker than those from the sneutrino mass limits given in 
equation ( p. 21 ). For the A" couplings, the indirect gluino contribution is far larger 
than the direct contribution, and we find that for a light spectrum with itlq and M1/2 
of order 100 GeV, the product of A^j-jAijsl^z) should be ~ 0.2. This bound scales 
very roughly as M1/2, but with some complicated dependence on the sbottom squark 
mixing as well as on the mass spectrum. 

4 — mixing 
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4.1 SUSY contributions 



The — sector has long been a probe of physics beyond the standard model, 
starting with the original motivation of the study of CP violation and recently also as 
a way of investigating flavour changing neutral currents in SUSY. These occur even 
in the absence of i?-Parity violation, since in the MSSM (and SM) quark flavour is 
not conserved. Hence we shall find that the RPV contribution is complementary to 
the already existing CKM contribution. We shall ask the question of how much it is 
possible to constrain RPV couplings through their impact on Aitik- Here we shall 
consider the relative sizes of three effects, the MSSM Am^ contribution and the direct 
and indirect RPV contributions. The first has already been investigated [[ill , |12|, as 



has the last [P, |], [T9|, but the question of whether the simple direct effects are 



in fact dominant has never been studied. 



The upper bound on the Kl — Ks mass difference has been measured |TJ] as 

ArriK = (3.491 ± 0.009) x lO'^^GeV (4.1) 
j{0 _ ^0 ]-Qixing is generated by the effective Lagrangian 

+ CRR[d,PRs'] [djPns^] + c'j,j,[d,Pns'] [d.Pns'] 
+ CLR[d,PLs'] [djPns^] + c'Ln[d^PLs'] [d.PRs'] 
+ dLL[d^i^PLs'][djYPLs']+dRR[da^PRs'][d,YPR^ (4.2) 

giving a contribution to the Kl — Ks mass difference 
Amft- = ^fx^riK {clr + "^Clr + ^^ll + ^dRR 



m 



^ {nid + msy ^^^^^ + + ^^^^ ~ ^LL + ^Crr - 4^) J (4.3) 

Here the terms involving both left and right handed fields in the effective lagrangian 
will give a larger contribution to Amx than those involving purely left or right handed 
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fields because of tlie factor of m^/ml. We shall not try to include QCD corrections to 
these formulae since these are unlikely to be more significant than the errors implicit 
in, for example, our choice of 0^3, m^, and SUSY parameters. 

Contributions to the parameters in the effective Lagrangian come from a variety 
of different classes of diagram, some of which are shown in Figure |[ The standard 
model contribution from Figure ^ is 

- 1287r2l'''="l ^ ^ ^ 

where for the second line we only consider the charm and up quark contributions. 
Apart from the standard model contributions, there are contributions from direct 



i?-Parity violating diagrams which have been calculated previously |T^, p2 
The only tree level diagram is that with the interchange of a sneutrino and two A' 
couplings shown in Figure Hb 



CLi.^' = -E^^#^ (4-5) 



m 



where our notation is as in the previous section and our earlier work ||T0 
The A' box diagrams such as Figure ^ give a contribution 

dLL^' = E 7^^^ife^'i2fcA'*i^A'i2m/4(m?^,m|,m^^,m^^) 

i,j,k,m 

+ E ■^^^'*iik^'j2k)^']im^'i2mh{m%,ml^,m^^^^ (4.6) 

i,j,k,m 

+ E -^^'iki>^%^'inii^\n,2h{rnl^,ml^ , m|, , m|„) 

i,j,k,m 

+ E ^^^'ifci^'ife2^W^'L2^4(m|,mJ,^,m2^,m2^) 

i,j,k,m 

+ E ^^^'ifci^'ife2^W^'L2^4(m^^,m^.,m2,,m|™) (4.7) 

i,j,k,m 
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i,j,k,m 



~ '^^^'iik^'j2k^Umi^'jrn2h{ml^,ml^,mlk^,m 



i,j,k,m 



while A" box diagrams such as Figure ^ give 



327r 



I r / 2 2 2 2 \ 



(4i 



(4.9) 



In addition we have diagrams such as those Figure ^ with one internal W boson 
and one internal squark line, where a helicity flip is needed on the internal fermion 
line, forcing it to be a top or charm. Hence we have 

\"W \ll \ll * T^* 7 f 2 nf2 2 2\ 

Clr = - 2^ -2^\ai\j-i2 Ki^Kj2mu,mu^J4{mi, M^, m„^, m ) 

• • 47r sm « 



/ \"w 

'^LR 



E 



47r sin^ 



A:3iA;32*^;i^.-2m„,m„^. J4(mf^, M^, m^^, (4.10) 



We neglect similar diagrams with A' vertices, since these are always dominated by 
the tree level contribution. 

Although they have not been considered in the context of i?-Parity, the diagrams 
involving mass insertions on squark lines have been analysed in the case where the off- 
diagonal mass terms are generated by the CKM matrix and by boundary conditions at 



the uniflcation scale [|TT], [1^, |18l . Most of these analyses include only gluino mediated 
box diagrams, such as Figure |^. 
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Am^ 7 



6^l^4«,m|^,M|,M|) 



-M?Jl'(m2- ,m2- M^) 

g 94V dj^' rf^' g5 gJ 



(4.13) 



(4.14) 

(4.15) 
(4.16) 
(4.17) 
(4.18) 



where we have assumed that mj^ ^ ""^s^- These results and those for the resulting 
contribution to Aitlk in equation ( |4.3| ) agree with those presented in Gabbianni et al 
1^ , where the disagreements between these results and others in references ITT], |12 



are discussed. 



The chargino mediated box diagrams give 



I'LL 



Ami 



a 



E I I ' I r^r« , < , ml , m|± ) (4.19) 



while from neutralino mediated box diagrams we have 
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(4.20) 
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(4.22) 
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(4.23) 



and from those with one gluino and one neutrahno : 



d 



RR 



^LR 



\Aau\''lK^i,mi,Mlo,Ml) 



(4.24) 



(4.25) 



(4.26) 



Here we use 



(AdRiAdLi + A; 

x/:(m|^,m|^,M|„,M|) 
+ {AdRiA*aLi + A*,j,,Aau)M^oM-,J':{ml,ml^, M|o, M?) 



■' cosBw \2 / 



(4.27) 



^dRj 



■' cos9"' ■' 
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and we note that 



A<<|, = (4.29) 

\2 oj/jo ^ -JO (5 \ ^ 



5mj 



\ dL da , 

^<,dA<nd. = ^1 %^^12 + ^^2lk2VQ,Q.^ + ^lVd,d.^ 

\ \ di^ dji, 

^ ^d ( ^d ^2 ^ I _d _2 



az, dn , 
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where the elhpsis stands for higher order contributions, the derivatives are assumed 
to act only of one of the arguments of the appropriate I'l, and we have only included 
terms proportional to Am? -, and Am? r in the dominant gluino contributions, not 
in the neutralino contributions where such terms also appear but are smaller. 

4.2 Analytical Discussion 

We begin by noting that the standard model term from equation (^4.4| ) gives a con- 
tribution to ArriK of around 2 x 10^^^ GeV, which given the large theoretical errors 
in the input parameters to the calculation is in reasonable agreement. Hence we 
shall here require that the new contributions from SUSY and RPV do not destroy 
this agreement by being larger than the experimental limit themselves. However, 
it is important to note that there are unknown relative signs between the various 
contributions, and hence that there can be cancellations. 



The direct contributions have been discussed in refs. p, ||, |I9|, The most 

stringent bound on couplings is that from the tree level diagram of Figure §3 which 
leads to a constraint [S, ^ 



X'^l2X'i2l{Mz) = X'ij2X'iki{Mz)Kj,Kk2 ~ 1.3 X 10-^ { ^q'^q^y )' (4-30) 

The box diagrams of Figures ^ and ^ and the competing diagram with an internal 
W line of Figure |^ lead to bounds on products of two couplings of order 10"^ to 
10^'^ 0, ^] for a very light spectrum. Ref. [|T^ was mostly concerned with the case 



where, in the weak basis, only one i?-Parity violating coupling is non-zero. This leads 
to bounds of order 0.1 on A'jifc and X'i2k for a very light spectrum, but scaling more 
weakly with masses. 

We now consider how large the expected contributions from the main indirect 
processes are, beginning with non-zero A", by deriving approximate bounds from 
each of the most important direct and indirect contributions in turn. The term which 
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we shall bound will be A"^ = A-3]^A^32, and for purposes of the discussion in this section 
we shall assume that all superpartners are degenerate at m ~ 3Mi/2- 

The main indirect contributions are the gluino mediated box diagrams, which 
arise because Am? j is non-zero. Using the various functions in Appendix C and 
Eqs. ( B.l ) and ( [4.26 ) we can find (ifj^ and hence derive an approximate bound of 



Although other indirect contributions exist, this is the simplest and is sometimes 
the largest. However, the various contributions from Am? r and Am? -r often in 



practice give comparably large effects, and so equation ( |4.31| ) should be treated with 
caution. 

The direct contributions are those from the A" box diagram and the A" — W 
diagram The former gives a bound of 



which may be rather inaccurate if the first index of the non-zero A" is 3. The second 
gives 

yi,,^'i,JMz) ^ <''<'5(^5^) (4.33) 

X'iMh) i 0.1 (^)^ (4.34) 

< f rh 



KMM.) ~ 0.2[^^) (4.35) 



A3'.A;.3(M.) ^ 0.1 {^)' (4.36) 
respectively, where the scaling is very approximate for the cases involving stop squarks. 



These numbers are in reasonable agreement with those in reference given the un- 
certainties in (for example) m^. 

We conclude that although the indirect contribution appears rather smaller, none 
of the contributions to Ami^ is obviously negligible in the region of interest, and so we 
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must turn to a numerical analysis, noting that we expect X"'^ (Mgut) ~ 10 ^ to give 
a result comparable with the experimental limit for masses of order a few hundred 
GeV. 

We now turn to the case of non-zero A'. Here the tree level diagram is inevitably 



completely dominant where it exists P, 22|, and even with only one non-zero coupling 



in the weak basis there can be measurable effects [0. Note that here it is possible 
to arrange the couplings so that the indirect contributions are zero but the direct 
contribution is not, for example if we arrange for non-zero A'113, A'223, A'212, and X'i22- 
However, such scenarios seem contrived given the need to avoid the extremely tight 
bounds from the tree level term while simultaneously having non-negligible values for 
four couplings, requiring remarkable cancellations in A'ji2A'j2i- 

The only scenario which we shall consider is thus non-zero X'n^ and A'j23- For this 
case the direct contribution lead to a constraint which is similar to that of Eq. (|4.31| ), 

X'^Mz) ~ 0.07 ( (4.37) 
^ ' VSOOGeV/ ^ ' 

where again this result is very approximate given the number of different contributions 

which can be relevant. There is no X' — W diagram by construction since it only exists 

when the tree level diagram exists, while the direct contribution from the box diagram 

then gives 

A-(M.)^8xlO-(^) (4.38) 

Since these results are in general only reliable as order of magnitude estimates, we 
use a numerical study to find out which contributions are most significant. 

4.3 Numerical Results 

The results for Attlk are fairly straightforward. We find that for each scenario which 
we have considered the indirect contributions are in fact completely dominated by 
the direct ones, with the gluino term being between one and two orders of magnitude 
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too small to compete with the direct box diagrams. We illustrate this in Figure ^ 
which shows a typical situation. Here we have set \2i^{Mgut) = ^223 (^gi/t) = 0.02, 
Ml/2 = 100 GeV, Aq = 0, tan/5 = 10, fi4 < 0, and show the relevant contributions. 

Apart from the fact that the direct contributions are dominant here, we note 
that the contribution from diagrams with one neutralino and one gluino line is com- 
parable in magnitude to that from gluinos alone. This clearly will be significant 
for models involving mass insertions from a GUT theory. In general we find that for 
flavour violation in the left handed sector the contributions from charginos and mixed 
gluino-neutralino diagrams are up to half that of the gluino, while for right handed 
flavour violation the chargino contribution is negligible but the mixed neutralino- 
gluino contribution is similar. 

In summary then, we find that in practice this is the only process which has 
been studied where the indirect contributions are generically overwhelmed by the 
direct one, and hence we can essentially simply use the results quoted above in equa- 
tions ( [4.32|) and (|4.38|) to bound the couplings, thus confirming the results of refer- 



ences 



191 , pi| Why this should be so is unclear, since the formulae for the 
contributions to AniK are so complicated, but it appears the values of the various 
four point functions are simply such as to suppress the indirect contributions, while 
for the three point functions relevant for /i — > 67 and b s'j the indirect contribu- 
tions are enhanced. We also find that different parts in the gluino contribution are 
of different signs and comparable magnitude, and so tend to partially cancel. 

Finally, we note that the qualitative nature of this result is such that the direct 
contribution will also dominate in 5° — and — mixing. 
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5 Conclusion 



We now conclude with a brief summary of our results. We have extended our analysis 
of the RGEs for i?-Parity violating supersymmetry to include the effects of the CKM 
matrix, and have studied two well-known flavour changing processes, including both 
the direct contributions, with i?-Parity violating couplings at the vertices of diagrams, 
and the indirect ones, where flavour violation arises through soft masses generated 
by the Renormalisation Group Equations. 

We first showed that constraints on A' can be derived by demanding that the 
sneutrino masses not be driven below their experimental limits, which gives a bound 
of around 0.3 on the low energy values of the A' couplings, with the particularly 
interesting feature that it is extremely insensitive to the sparticle spectrum and does 
not disappear in the limit where the masses go to infinity. Given that the sparticle 
spectrum is now known to be heavy (with many limits well above 100 GeV if we 
assume unification) this means that this bound is one of the tightest existing in such 
a scenario. 

For the case of 6 — > 57 we find that the indirect contributions caused by the flavour 
violation are dominant, and enhance the amplitude contribution by up to an order of 
magnitude, allowing new constraints on both A" and A'. However, these are rather 
weak, and the A' constraints are weaker than those from sneutrino masses derived 
earlier in this paper. 

With regard to — mixing, we have included all the indirect contributions 
including those from gluino, neutralino, chargino, and mixed neutralino-gluino dia- 
grams. We find that here the chargino and mixed diagrams are sometimes comparable 
in size to those from gluinos, with consequences for the study of flavour changing from 
GUT theories. However, in the context of i?-Parity violation, we find that the indirect 
contributions are consistently around an order of magnitude or more less than the 
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direct ones, so that for this process and for the similar ones of — and — 



mixing we cannot improve on results in the literature [19], . 
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A Renormalisation Group Equations 

We now explain how to derive the RGEs which we shall use. Full RGEs for our model 



are given in reference but there are two slight complications which arise when 
considering quark flavour violation which do not occur in the lepton flavour violating 
processes considered there. The first of these is that the MSSM (and the SM) contain 
QFV explicitly through the CKM matrices and through the MSSM soft trilinears. 
The second is that the RGEs of are presented in a weak basis while we work in 
a fermion mass eigenstate basis. 

The inclusion of the CKM matrix is straightforward, since defining 

= 5,,h'^ (A.l) 

with hf one of the eigenvalues of the Yukawa matrix /if^-, and similarly for A**, A*^, we 
can simply put 

X\^- = -hf- = -K*A'^ = -A'^ - kA'^ (A.2) 

where K is the CKM matrix, which we shall henceforth consider to be real, defining 
K = 1 + k so that k is anti-symmetric to first order in k. A full description of our 



other notation is contained in reference [10|. This leads to off-diagonal terms in the 
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running of the soft squark masses 



''i 



^ 4fc,,/i:^/i,^K^ - + 4r/f,<. (A.3) 

which are in addition to the i?-Parity violating contributions. Here rj^ is given at the 
unification scale by Vj^ = h'^Ao = K*A'^Aq, but in general need not run in such a way 
as to preserve this relationship, and we have assumed that the diagonal terms in the 
mass matrix := m^^ are very much larger than the off-diagonal ones. Including 
the down and lepton soft trilinears in the equations given in ||10| can be done by 
treating Hi noted there. 

Including the full matrix of rj^j is straightforward. We give here the full RGE 



-1-4/)" X" C" -I- 9r)" X" X" -I- r)" X' X' -i-97)"X' C 

'^"'i ^imn^jmn ^' lik^kran^jmn 'Ikj^ mkn^ mm ~r ^it'j ^ mjn^ m 



/13 2 o 2 16 2 
+25,,/ir {^^9\Mi + ?>glM^ + y^?|^3) (A.4) 

where we have included the down and lepton Yukawa couplings and trilinear soft 
terms only implicitly [i.e. through the A, A', C, and C'\ We also must then include 
a term in the RGE for C[^^ of 2ri^ihfX'uk + hfClj^ and one in the RGE for C^'j^ of 
Arjf^h^X'ljf, + 2h'!^'^C"ji^. Inclusion of the remaining Yukawa couplings is trivial. 

The second and rather more complicated issue is that of the difference between 
the weak and mass bases. We have three Yukawa matrices /i", h'^, h^, and it is 
conventional to choose a basis such that 

/i« = A" h'^ = KK'^ h^ = A^ (A.5) 

where K is now being taken real. Such a choice is also a weak eigenstate basis. 
However, these relations are not stable under the RGEs with either a non-trivial 
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CKM matrix or i?-Parity violating couplings. Hence we will find in general that all of 
the Yukawa matrices will have non-zero off-diagonal elements, greatly complicating 
our analysis. The simplest way of dealing with this is to define a basis where the field 
rotations to impose equation ( |A.5| ) are performed in a scale-dependent way. Then 



equation (|A.5|) will always be satisfied, but the RGEs will differ slightly. 
If we define 

A" = t//i"Vf (A.6) 

where U, V are orthogonal matrices with explicit dependence on the renormalisation 
scale, then 

A" = f7/i"Vf + t//i"Vf + t//i"Vf (A.7) 
where the dots indicate IGn'^d/dt and defining u = UU^ and v"^ = V^Vj we find that 

A" = mA" - AV + 7'3A" + A" (7" + 7^2) (A.8) 

where the anomalous dimension matrices 7 are solved in the basis where equa- 
tion (|A.5|) is satisfied. Hence we can solve for u and in terms of the Yukawa 
couplings and anomalous dimensions. The situation for the down-type Yukawas is 
slightly more complicated. We choose to define 

A^ = K^iWh'^V^) (A.9) 

where K is the CKM matrix and W and are unitary as before. Hence we may 
calculate w and defined hj w = WW^ and v'^ = VdVj'. For the leptonic case we 
define and v"^ exactly as for up-type quarks. 

Imposing that the A remain diagonal we find 

_ _ Q ( hf + hf \ _ „ ( 2h^h- \ 
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hf - hf 

Here we use w = K^wK and 7 = K^'yK. 

Explicitly expanding taking the masses to be very strongly ordered with hi <C 
^2 ^ ^3, we find that for i > j 




V 



mn ^2 mn 

9_ii. 7,°'^ _|_ 9X" X" -|-9^V^X' X' 

mn mn 



h 7)"2 1 \/ \' _\A^\^(\' X' -I- X" X" 

jj rojj/tj -r ^ /\ min^ mjn "r ^ r mni^ mnj T ^mni^mnj i 



o_lu 7,"2 _|_ 9 v^/A^/ \/ . \/f \/f ^_|_9^V^x' X' 

l.d V i mni^ mnj ~T ^m.ni^mnj ) ~^ r^d / ^ min^ mjn 

i mn i mn 

^^(^^imn^jmn ~l~ 3A jmn-^ jmn) ~l~ 4 ^ ^ ■^mni'^mn 



/ A ^^^mn'^imn i '-'^^ imn^^ jmn) i / ^ ^^mm'^mni 

mn mn 

2 ^ ^ ^mni^mnj ~l~ '^ITe 5Z "^J™" ~^ "^"^ imn-^ jmn) 
mn i mn 



-%(/^f + /^f)+4^E^:mnA;n 



"i mn 



~^Td E^"^ mm A mnj + ^mni^mnj) (^-H) 
"'i mn 

Hence we can simply use the RGEs in reference |TD[ , but adding the explicit CKM 
contributions and including the conversion of basis by, for example 

, 9 mass basis , o weak basis 

This gives the fermion mass basis RGEs for the off-diagonal terms in the soft 
masses 

16vr'— 1^ = E(A'--A'™."('^L+"^Q.+2< + 2mi^) + 2C^„C;,„) 
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+2kij{hf - hf)m^ 



Hi 



+2nt,ril + 2niri% + u,^[m\ - m|J (A.13) 

mn 

+4« + ^.';«-<,) (A.14) 

mn 

+AhXh^{ml^ - mi;) + ^ri-ni^ + v%{m\ - mD (A.15) 

mn 

+ T.{^X'^mnyjmn{ml + 171% + 2ml + ^KJ + ^CLnC'jmn) 
mn 

+HkV-k + ut,{m%-ml) (A. 16) 

1677^—^ = Y,{Xrnni>^mnj{ml. + ml. + 4m|^) + 2CmniCmnj) 

mn 

+4«,- + (A.17) 

It is clear from these equations that when we consider the first and second generations, 
where the soft masses are nearly degenerate, the difference in basis is totally irrelevant, 
but for cases involving the third generation which is typically rather lighter the effects 
may be non-negligible. 

The RGE for the CKM matrix is then 

dk 

167r^ = Uij-Wij (A. 18) 

B Approximate Solutions of RGEs 

For our discussion of the indirect contributions to the various processes, it will now 
be helpful to present some approximate formulae for the various ofi^-diagonal mass 
insertions. For convenience we will present here formulae for all the elements, not 
merely those involved in 6 ^ S7. If we begin by switching off the CKM matrix we 
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find 

ml.Q. {Mz) ^ >^'minX'mjn{MGUT) (5mg + ISM^^ - 5^oMi/2 + Al) 

mn 

mldj (Mz) = X'mmX'mnj{MGUT){lOml + SOM^^ - lOAoMi/2 + 2Al) 

mn 

- E >^'Li^'LjiMGUT){lOml + 50M,% - 20^Mv2 + 5^^) 

mn 

mlu,iMz) = -Y.X'L^X';^^{MGUTmml + 50M'y,-20AoM,/2 + 5Al) 

mn 

(B.l) 

while for the Am|^ and Amj^j^ insertions we find 

Vu^^l = m,Y.^'mlnX'm2niMGUT)(^Mi/2-l.5Ao) 

mn 

rj^^ui = m,^A'„„iA'^„2(MGc/T)(9Mi/2-4.5A) 

mn 

+ rrisY: \'LAAMGUT)my2 - 3.5Ao) (B.2) 

mn 

where is the running strange mass at Mz- Since all these formulae depend quite 
strongly on a^, they should not be trusted to better than about a factor of two. Note 
that the Am|,^ and Am^^ insertions are very much smaller than the others because 
they are proportional to a soft mass times a strange quark mass rather than two soft 
masses. 

Including the CKM matrix is straightforward, and leads extra contributions to 
those given above of 

^Q.Q.iMz) = — (3m^ + 10Mi%-4AoMi/2 + ^g) 

COS kJ 

^QiQsW = ^(2m2 + 8Mi%-3^Mi/2 + 0.8^g) 

COS lj 



^Q^Q^i^z) = — (2m^ + 8Mi^/2-3AoMv2 + 0.8A^) 

COS l) 

?7iVi = m,(-0.8Mi/2 + 0.2^) 

« r^l^u^ (B.3) 

although here the errors are larger, owing to the dependence on poorly known Yukawa 
couplings. Although we have only shown the contributions to ttIq.q., in fact there 
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are also contributions to ^^2^3 ^U2U3 particular which are much smaller but 
non-zero. 



C Definition of Functions 



The expressions for the amplitudes in 6 ^ S7 and K K mixing presented earlier 



employ a number of functions, many of which are contained in ||T0[, and the remainder 
of which are given explicitly here. Functions occurring in 6 — 57 are 



J n 



(2), 



X 



J a 



f(3) 



J a 



(4), 



X 



f(4: 



(5), 



X 



f{5 
J a 



(6), 



(6), 



X 



X 



QuFiix) + F2{x) 
Fi{x) 

QuFsix) + F^ix) 
Fsix) 

F,{x) + Q^F^ix) 
F2{x) 

Fsix) + QuF^ix) 
Fi{x) 

-QuFiix) + QdF2{x) 
^F,{x) - ^F2{x) 
-QdFi{x)+Q^F2{x) 



(C.l) 



In addition, for the ^uik calculation the functions 14 and J4 used in the text are 
given by 



Til 1 1 1\ 

iii\m^, mg, rrig, vrii^) 

7/2 2 2 2\ 

where we have defined 



dxdydz — ^ 

-D4 



(C.2) 



D4,(ml, ml, ml, ml) = [xm\ + (1 — x)m^z + [ym\ + (1 — y)m^^{l — z) (C.3) 
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Note that I4 and J4 are totally symmetric about interchange of arguments, and can 
be simplified in certain limits as follows. 



74(m^m^M^M2) 
himl,ml,M^,0) 

h{mj,ml,M^,M^) 

/4(m2,M2,0,0) 
74(m^m^M^0) 

I'l{m^,m^,Ml,Ml) 
I'l{m\,ml,M\M^) 



J^{mi,mi,M\M') 

74(0, mi, 7712, M^) 

J'l{rn^,m^,M'^,M'^) 
J'l{mlmlM'',M'') 



1 
1 



1 

1 

2m^ 



(1 — a;i)a;2 lna;2 — (1 ~ X2)xi \nxi 
2(xi - l)(xi - X2)(a;2 - 1) 
1 



2(1-Xi)(l-X2) 

1 / x|ln(a:2) xl\n{xi) 



X2-X1 V2(l-X2)2 2(l-xi)2 



In a; 



2(,x - 1) 
l-y + y\n{y) 



m4(Mf - M|) 
1 



yiFiivi) -Z/2-^i(z/2) 



M6 



Xi + ./•2 + .rf + .1-2 - G:r 1:1-2 + ■ri.r2 + .fixl 
2(a;i - l)2(a;2 - l)2(a;i - 2:2)2 



(,Ti,T2 — xf) ln(a;i) {xiX2 — x^) ln(a;2) 



{Xi - l)^{Xi - X2)^ {X2 - l)^{Xi - X2)^ 



1 



1 

2 



-1 



1 

M8 



(1 - Xi){l - X2) 

1 / In(xi) X2 ln(a;2) \ 
1X2 - xi \(1 - (l-a;2)V 
(2:2 - 1) ln(a;i) - {xi - 1) ln(a;2) 
(x2 - l)(x2 -xi)(a;i - 1) 



m° 



-2(1 — — a;2 + ^1 + X2 — X1X2) 

{xi-mx2-iy{x,-x2y 

{xi + X2 + X1X2 - Sxj) \n{xi) 

(,Ti - X2)^{Xi - 1)3 

(xi + .)-2 + .ri.r2 - 3.r^)lii(x2) 

{Xl - X2)^{X2 - If 

y2ln(y2) yiMvi) 



(y2 - mvi - y2) (yi - myi - ^2) 

ll + 7(yi + y2) 



6{yi - 1)3(1/2 - ir 
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(C.4) 
(C.5) 



(C.6) 

(C.7) 

(C.8) 
(C.9) 
(C.IO) 



(C.ll) 



(C.12) 

(C.13) 
(C.14) 



(C.15) 



(C.16) 



-2(2/1 + 2/2 - %i2/2) - 5(2/1 + 2/2)2/12/2 + 2/^2/2) 



Here x = m^/M'^, Xi = mf/M'^, y = M'^/m^ and 2/i = Mf /m?^ and 

Til f 222 2n r/222 2\ 

J" {ml, ml, ml, ml) = -^-^^^Ji{m\,ml,ml,ml) (C.17) 

We give /e and /g defined in reference [|rT| for convenience when comparing our results 
with other authors. 
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D Figure Captions 



Figure [l| 

Absolute values for various contributions to the b ^ sj amplitude Am. Parameters 
are X'^iiMcur) = A'131 (Mgc/t) = 0.05, mo = 100 GeV, Aq = 0, tan/3 = 10, /i4 < 
0. Contributions from Higgs (dotted), SM (horizontal dotted), neutralino (double 
dashes), chargino (dashes for the MSSM part with triple dashes for the remainder), 
gluino (quadruple dashes), direct RPV diagrams (dot dashed) and total (solid) are 
shown. The experimental upper and lower limits are shown as horizontal solid lines. 

Figure |^ 

Absolute values for various contributions to the 6 ^ 57 amplitude A^l. Parameters 
are X'Ii2{Mgut) = XiisiMcuT) = 0.1, mo = 100 GeV, Aq = 0, tan/5 = 10, /i4 < 
0. Contributions from neutralino (double dashes), gluino (quadruple dashes), direct 
RPV diagrams (dot dashed) and total (solid) are shown. The experimental upper 
limit shown varies with M1/2 because the MSSM part of Alr does. 

Figure |^ 

Some of the diagrams contributing to — mixing. 
Figure |^ 

Absolute values of various contributions to ArriK- Parameters are X213{Mgut) = 
X223{Mgut) = 0.02, Ml/2 = 100 GeV, Ao = 0, tan/3 = 10, /i4 < 0. Contributions 
from neutralino (double dashes), gluino (quadruple dashes), mixed gluino-neutralino 
(dashed), direct RPV diagrams with and without W lines (dotted and dot dashed 
respectively) and total (solid) are shown, together with the experimental upper limit 
(horizontal solid line). 
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